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We analyse H(z) parameter data with some conditions by using Lagrange mean value theorem 
in Calculus. We find that: (1) there exists at least one decelerated phase at 1 ø confidence level 
in the redshift range (0.38, 0.59); (2) the equation of motion of dark energy may be less than —1 
at 1 ø confidence level at some redshifts in the redshift range (1.3, 1.53); (3) there exists at least 
one accelerated phase at 1 ø confidence level in the redshift range (1.037, 1.944). These results may 
provide possible evidences for physics beyond ACDM. 


PACS numbers: 


I. INTRODUCTION 


A great number of independent cosmological observations, such as supernova Ia (SNIa) at high redshift [1, 2], 
large-scale structure [3], and the cosmic microwave background anisotropy [4, 5], have confirmed that the Universe 
is experiencing an accelerated expansion. In order to explain this phenomenon, an unknown energy component 
(dubbed as dark energy) usually have to be introduced in the framework of general relativity. The simplest and 
most theoretically sound scenario of dark energy is the vacuum energy (ACDM) with a constant equation of state 
(EoS) wx = px/px = —1 where wx denotes the EoS of dark energy. This model is consistent with most of the 
current astronomical observations, but suffers from the cosmological constant problem [6] and age problem [7] as well. 
Recently, Hubble tension may also provide evidences for physics beyond ACDM [8]. 

The general approach to studying dark energy is to assume either a theoretical model or an EoS, and then use 
observational data to limit relevant parameters, see for example, for spatially-flat ACDM the Hubble constant and 
the matter density parameter are constrained as: Ho = (67.4 + 0.5) km s~'Mpc~!, Qm = 0.315 + 0.007, respectively; 
while for EoS (wy = wo + 722) parameterized model, the related parameters are limited as: Ho = (68.31 + 0.82) 
km s~!Mpc™!, wo = —0.957 + 0.080, and wa = —0.2979'3% [5]. Statistical methods, such as the maximum likelihood 
(7, 9-12], are generally used to analyze the observational data to fit the parameters. These statistical method yields 
the best statistical results, but it is easy to eliminate some interesting (possibly important) data. Here we propose 
a model-independent method by using the Lagrange mean value theorem to analyze H(z) parameter data. We find 
that the EoS of dark energy may be less than —1 at some redshifts and the accelerated phase may occur earlier than 
we previously thought. 

The paper is organized as follows. In the next Section, we will present H(z) parameter data and derive the equations 
needed to analyze these data. In Sec. HI, We will provide the data and results obtained from the analysis. Finally, 
we will briefly summarize and discuss our results in section IV. 


II. THEORETICAL METHOD AND H(z) PARAMETER DATA 


In this Section, we will present 63 H(z) parameter data obtained recently, then introduce the Lagrange mean value 
theorem and combine it with Friedmann equations to derive equations needed to analyze H(z) parameter data. 
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A. Theoretical method 


Assuming a Friedmann-Robertson-Walker-Lemaitre (FRWL) spacetime 
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ds? = —dt? + a? (t) | —_—~ 
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where a(t) is the scale factor, K denotes the curvature of the space with K = +1, 0, and —1 corresponding to a 
closed, flat and open universe, respectively. We use the unit c = 1 here. According to the Planck 2018 results, 
the spacetime is spatially flat: Qko = 0.001 + 0.002 [5]. So we consider a spatially flat FRWL spacetime here, the 
Friedmann equations take the form 
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H? = — 2 
3 Ps (2) 
a 4TG 
A +3), (3) 
or equivalently 
H = —4nG(p + p), (4) 


where the H = å/a is the Hubble parameter with the dot denoting the derivative with respect to the cosmic time t. 
The total energy density p and pressure p contain contributions coming from the radiation, nonrelativistic matter, 


and other components. Because dz = —(1 + z)Hdt, we have 
. dH 
Combining Eqs. (4) and (5), yields 
dH 4nG AnGip(1+ w) 
= } = 5 6 
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where w, is the total EoS. From this equation, we can judge whether the total EoS is greater than, equal to, or less than 
—1: see for example, if dH/dz < 0, we have wx < ws < —1 because of the positive of H and p. In an era dominated 
by dark energy, we can also determine with Eq. (6) wether the EoS of dark energy is equal to —1: if dH/dz = 0, 
then one have wy ~ w, = —1. If dH/dz < 0, we know the Universe is experiencing an accelerated expansion. But 
if dH/dz > 0, we can’t judge whether the Universe speeds up. At this point, we need another important physical 
quantity, the deceleration parameter, which is defined as 


1 dH 
=-14+(1 ——. T 
q=-14+ (+e) (7) 
Now a question naturally rise: if we have some H(z) parameter data, how can we use them to directly determine 
dH/dz or q? Think of Lagrange mean value theorem in Calculus, which states: for a continuous and differentiable 
function f(x), there exists x, < £12 < £2 satisfying 
f(z2) — f(x) 


df 

dx = = ` (8) 
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Applying this theorem to Hubble parameter which we assume is a continuous and differentiable function of z, and 

taking function H(z) as f(x) in (8), we have 

dH H(z) — H(z;) 

H (ag) = A pay = a (9) 
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where zj < zij < 2. If zi—zj « 1, H’(z;) will be large in general, which will make relevant results less credible. Since 

zij and H(z,;) in Eq. (7) are unknown, we take approximatively: zi; = (z; + 2;)/2 and H(z;) = [H (zi) + H(z;)]/2, 

which can be called as mid-value approximate method. Then we have 


(2 + 24+ zj) H' (zij) 
H(z) + H(z;) 


q(zij) > -14 (10) 


index z H(z)[km s~'Mpc~'] om Reference Method index z H(z) on Reference Method 
zı 0 74.03 1.42 [8] SN Ia/Cepheid|z33 0.51 90.4 19 [I3 Clustering 
Z2 0.07 69 19.6 |14 DA z3a 0.52 94.35 2.65 [i5 Clustering 
23 0.1 69 12 [16 DA z35 0.56 93.33 2.32 [15 Clustering 
Z4 0.12 68.6 26.2 [14 DA zs 0.57 929 7.8 17 Clustering 
25 0.17 83 8 16 DA z37 0.59 98.48 3.19 [15 Clustering 
26 0.1797 75 4 18 DA z3g 0.5929 104 13 18 DA 

Z7 0.1993 75 5 18 DA z39 0.6 87.9 6.1 19 Clustering 
Z8 0.2 72.9 29.6 [14 DA zoo 0.61 97.3 2.1 13 Clustering 
Z9 0.24 79.69 2.65 [20 Clustering ZA 0.64 98.82 2.99 [15 Clustering 
zio 0.27 77 14 [16 DA za2 0.6797 92 8 18 DA 

zıı 0.28 88.8 36.6 [14 DA zag 0.73 973 7 19 Clustering 
Z12 0.3 81.7 6.22 [21 Clustering ZA 0.75 98.8 33.6 [22 Clustering 
zı3 0.31 78.17 4.74 [15 Clustering za5 0.7812 105 12 18 DA 

Z14 0.34 83.8 3.66 [20 Clustering Z46 0.8754 125 17 18 DA 

zis 0.385 82.7 8.4 [23 Clustering za7 0.88 90 40 16 DA 

zig 0.3519 83 14 [18 DA zas 0.9 117 23 16 DA 

Z17 0.36 79.93 3.39 [15 Clustering ZA9 0.978 113.72 14.63 |24 Clustering 
zig 0.38 81.5 1.9 [13 Clustering z50 1.037 154 20 18 DA 

zig 0.3802 83 13.5 [25 DA z51 1.23 131.44 12.42 [24 Clustering 
z20 0.40 882.04 2.03 [15 DA z5s2 13 168 17 16 DA 

za1 0.4 95 17 [16 DA z53 1.363 160 33.6 [26 DA 

z22 0.4004 77 10.2 [25 DA z5a 1.43 177 18 16 DA 

z23 0.4247 87.1 11.2 [25 DA Z55 1.526 148.11 12.71 [24 Clustering 
Z24 0.4293 91.8 5.3 [25 DA zs 1.53 140 14 16 DA 

zon 0.43 86.45 3.68 [20 Clustering z57 1.75 202 40 16 DA 

z96 0.44 82.6 7.8 [19 Clustering zss 1.944 172.63 14.79 [24 Clustering 
zor 0.44 84.81 1.83 [15 Clustering z59 1.965 186.5 50.4 [26 DA 

zag 0.4497 92.8 12.9 [25 DA zeo 2.3 224 8 27 Clustering 
Z29 0.47 89 34 28 DA 261 2.33 224 8 29 Clustering 
zzo 0.4783 80.9 9 25 DA 262 2.34 222 7 30 Clustering 
231 0.48 87.79 2.03 [15 DA 263 2.36 226 8 31 Clustering 
z32 0.48 97 62 [16 DA 


TABLE I: Hubble parameter compilation from cosmic chronometers (DA) or from the radial BAO surveys (clustering). 


If z; — zj is large, in general, Eq. (10) will be not valid. Taking the uncertainty on the values of H(z) into account, 
the uncertainties associated to the H’(z) and q(z) are given by, respectively 


and 


With oy and og, we can determine whether the results are credible at 1 ø confidence level. 
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index H' (z) jog q(z) |oq (index H' (z) low q(z) |oq 
zeı € (0,0.1797) |5.398 |23.62 |-0.921)/0.348) za714 € (0.34,0.44) {10.1 40.92 |-0.833|0.679 
z71 € (0,0.1993) |4.867 |26.08  |-0.928]0.387) z359 € (0.24, 0.56) 42.625 |11.006 |-0.31 /0.192 
zo1 € (0, 0.24) 23.583 |12.527 |-0.656]0.189] z4o9 € (0.24, 0.61) 47.595 |9.138  |-0.233)/0.162 
z131 € (0,0.31) 13.355 |15.962 |-0.797/0.249]za19 E€ (0.24, 0.64) 47.825 |9.988  |-0.228]0.178 
z1a1 E (0, 0.34) 28.745 |11.547 |-0.574| 0.182] z3418 € (0.38,0.52) [91.786 |23.291 /0.514 |0.412 
z171 € (0, 0.36) 16.389 |10.209 |-0.749|0.162} 23420 € (0.40,0.52) |102.583 |27.818 |0.698 |0.493 
z181 € (0, 0.38) 19.658 |6.242 |-0.699|0.1 | z3718 € (0.38,0.59) /80.857 |17.68 0.334 |0.319 
z201 € (0, 0.4) 20.025 |6.193  |-0.692|0.1 | z3431 € (0.48,0.52) [164.00 |83.454 /1.701 |1.424 
z251 € (0.0,0.43) |28.884 |9.173  |-0.562]0.15 | 23727 € (0.44,0.59) [91.133 |24.518 |0.507 |0.436 
zo71 € (0.0,0.44) [24.5 5.264  |-0.624| 0.086] z3731 E€ (0.48,0.59) |97.182 |34.374 |0.602 |0.599 
z311 € (0.0,0.48) |28.667 |5.161 |-0.561]0.086) 24324 € (0.4293, 0.73)|18.291 18.959 |-0.694| 0.326 
z139 © (0.24,0.31) |-21.714|77.578 |-1.351] 1.241) 25150 € (1.037, 1.23) |-116.891] 121.983 |-2.747| 1.679 
zi7g E (0.24, 0.36) |2.00 35.857 |-0.967|0.585 | zs550 € (1.037, 1.526) |-12.045 |48.46  |-1.182]0.718 
zig9 € (0.24,0.38) |12.929 |23.291 |-0.79 |0.383) 25650 € (1.037,1.53) |-28.398 |49.519 |-1.441|0.733 
z209 € (0.24,0.40) |14.688 |20.864 |-0.76 |0.346) 25552 € (1.3,1.526) |-88.009 |93.92 |-2.344|1.344 

€ (0.24, 0.43) 135.579 |23.868 |-0.428]0.399] zs652 € (1.3, 1.53) -121.739]95.75 — |-2.909) 1.365 
z279 © (0.24,0.44) 25.60 [16.102 |-0.583]0.271) 25554 € (1.43, 1.526) |-300.938] 229.532 |-5.588| 3.188 
z319 E€ (0.24,0.48) |33.75 | 13.9091]-0.452) 0.237] 25654 € (1.43,1.53) |-370 228.035 |-6.789 | 3.152 
z339 © (0.24,0.51) |39.667 |12.077 |-0.359] 0.208) z5g54 € (1.43, 1.944) |-8.502 45.324 |-1.131|0.688 
z2014 € (0.34, 0.40) |-29.333]69.755 |-1.485 1.14 


TABLE II: H'(z) and q(z) data obtained from H(z) parameter data. 


B. H(z) parameter data 


The data set we use consists of 1 H(z) measurement from SNIa observation, 34 H(z) measurements obtained by 
calculating the differential ages of galaxies, which is called cosmic chronometer, and 28 H(z) measurements inferred 
from the baryon acoustic oscillation (BAO) peak in the galaxy power spectrum, as listed in Table I. In three cases, 
the datasets are given with their lo confidence interval. 


II. APPLICATIONS 


In this Section, we apply Eqs. (9), (10), (11), and (12) to investigate the evolution of the Universe with the 
observational Hubble parameter data. In order to avoid significant errors, we have considered the following limitations 
in the process of analyzing the H(z) data: 0.1 S z; — z; $0.5, og < 5 if H < 100, and oy < 20 if H > 100. The 
data for H’(z), q(z) at 1 o confidence level are listed in Table I, from which we can conclude: 

(a) At redshifts 61, 271, 291, 2131, 2141, 2171, 2181, 2201, 2251, 2271, 2311, 2139, 2179, 2189, 2209, 2259, 2279, 2319, 2339, 
29014, 227145 2359) 2409; 2419) 24324, 25150, 25550) 25650) 25552, 25652 25554, 25654, aNd Z5g54, the Universe experiences an 
accelerated expansion at 1 o confidence level. 

(b) At redshifts 23418, 23420, 23718, 23431, 23727, and 23731, the Universe experiences an decelerated expansion at 1 
o confidence level. 

(c) At redshifts 2139; 220145 25150) 75550) ~5650) 255525 and 25854; since H'(z) < 0, implying Wx < Wg < —1, but not 
at at 1 ø confidence level. However, we can infer that wy < wą < —1 at redshifts 25652, 25554, and z5654 at at 1 o 
confidence level. 

According to the Planck 2018 results [5], the matter density parameter for the spatially-flat ACDM was constrained 
as: Qm = 0.315 + 0.007, implying that the phase transition from deceleration to acceleration of the Universe occurs 
at the redshift z ~ 0.632. Result (b), however, shows that there exists at least one decelerated phase in the redshift 
range (0.38,0.59) at 1 ø confidence level. In addition, results (a) and (c) also indicate that the EoS of dark energy 
may be less than —1 at some redshifts and the accelerated phase may occur earlier than we previously thought. 


IV. CONCLUSIONS AND DISCUSSIONS 


Using the Lagrange mean value theorem in Calculus, we have analysed H(z) parameter data with some conditions 
and find that: (1) there exists at least one decelerated phase at 1 ø confidence level in the redshift range (0.38, 0.59); 
(2) the EoS of dark energy may be less than —1 at 1 ø confidence level at some redshifts in the redshift range 
(1.3, 1.53); (3) there exist at least one accelerated phase at 1 ø confidence level in the redshift range (1.037, 1.944). 
The q(z) data we obtained can be used to investigate cosmological models. The results may provide clues to address 
Hubble tension and possible evidences for physics beyond ACDM. 
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